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ABSTRACT 

The a turbulent viscosity formalism for accretion discs must be interpreted as a mean 
field theory. The extent to which the disc scale exceeds that of the turbulence determines 
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\Q ' the precision of the predicted luminosity L u . The assumption of turbulence and use of a 
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implies: (1) Field line stretching generates a magnetic pressure > a /6 of the total pressure 
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generally, and a 1 to 1 relation between a and the pressure ratio when shearing instabilities 
Qh 1 dominate the viscosity. (2) Large eddy sizes and speeds in typical advection dominated 



accretion flows (ADAFs) lead to a lower precision in L v than for thin discs of a given 
total observation duration and central mass. The allowed variability (relative precision) at 
a particular frequency increases (decreases) with the size of the contributing region. For 
X-ray binary type ADAFs, the allowed variability is ~ 5% at R < 1000 Schwarzchild radii 
for averages over > lOOOsec. But for large galactic nuclei like NGC 4258 and M87, the 
relative precision error can approach 50 — 100% even at R < lOORs for currently available 
observation durations. More data are then required to compare with ADAF predictions. 

Key Words: accretion discs; galaxies: nuclei, active; Galaxy: centre; turbulence; binaries: 
general 



1. Introduction 

Accretion discs are a widely accepted paradigm (e.g Pringle 1981; Papaloizou & Lin 
1995) to explain a variety of features in high energy astrophysical sources such as active 
galactic nuclei (AGN), X-ray binary systems, cataclysmic variables (CVs), and dwarf no- 
vae. As accreting gas orbits a central massive source, internal energy dissipation drains the 
rotational energy, allowing material to move in and angular momentum out. The dissipa- 
tion sustains steady accretion and some fraction of the dissipated energy accounts for the 
observed luminosity. Micro-physical viscosities are too small to explain observed luminosi- 
ties so an enhanced transport mechanism, likely involving turbulence, is essential. Since 
astrophysical discs are surely magnetized to some non-zero level, the "Balbus-Hawley" 
shearing instability (e.g. Balbus & Hawley 1991; c.f. Balbus & Hawley 1998), which 
produces self-sustaining turbulence, is a natural and likely ubiquitous driver of angular 
momentum transport for at least thin discs and possibly thick discs as well. Significant 
dissipation may also occur above a thin disc, in a corona (e.g Haardt & Maraschi 1993; 
Field & Rogers 1993; DiMatteo et al. 1997). For low enough accretion rates, the dissipated 
energy may be primarily advected rather than radiated (Ichimaru 1977; Rees et al., 1982; 
Narayan & Yi 1995ab) forming an advection dominated accretion flow (ADAF) thick disc. 
ADAFs have been effective in modeling quiescent phases of accretion in a variety of X-ray 
binaries and galactic nuclei (see Narayan et al. 1998a, for a review). 

While non-linear instabilities in thin discs have been extensively simulated locally (e.g. 
Brandenburg et al 1995; Stone et al., 1996 Balbus et al., 1996) a useful approach to global 
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disc models has been to swipe the details of the stress tensor into a turbulent viscosity of 
the form (Shakura 1973; Shakura & Sunyaev 1973) 

Vfb = ac s H ~ vtbkb, (1) 

where H is the disc height, c s is the sound speed, l t b is the dominant turbulent eddy 
scale, v t b is the eddy speed at that scale, and a < 1 is taken to be a constant. Use of 
this formalism requires a mean field theory. Viscous coupling of differentially rotating 
fluid elements is a local paradigm, so assumptions of azimuthal symmetry and steady 
radial inflow (e.g. Pringle 1981; Narayan & Yi 1995ab) require turbulent motions to be 
statistically smoothed over the time and/or spatial scales on which mean quantities vary. 

The required mean field approach is similar to that employed in mean-field magnetic 
dynamo theory (Parker 1979), where the induction equation for the magnetic field is aver- 
aged and solved. In the kinematic (and therefore incomplete) dynamo theory, the velocity 
is imposed and the momentum equation is ignored. For the simplest global a accretion disc 
approach, the focus is reversed (but also incomplete); the momentum, energy and continu- 
ity equations are solved, with the inclusion of the magnetic field as a pressure rather than 
employing the magnetic induction equation (e.g. Narayan & Yi 1995ab). However, the 
usual disc equations are not usually derived from the formal averaging approach. Balbus 
& Hawley (1998) have addressed some of these points, but the conditions for which the 
standard equations result when H ~ R are not commonly studied. 

Despite being incomplete, the a formalism provides a useful framework for thin and 
ADAF discs. The absence of much scale separation between H and R has important 
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consequences. Here I show that l t b < H and field line stretching place a non-trivial lower 
limit on the magnetic energy for large a. A 1-to-l relation between a and (3 results when 
the viscosity is due to shearing instabilities as shown in section 7. I also estimate the 
precision of the a turbulent disc formalism and interpret this for thin and thick ADAF 
discs. For ADAFs the precision is lower, and thus the allowed variability higher than for 
thin discs with a given central mass, and fixed observation duration. A discussion of the 
implications for stellar vs. galactic nuclei presumed ADAF systems is given. 

Balbus & Gammie (1994) present a nice discussion of fluctuations in a thin disc and a 
relation to luminosity, but a different approach and different questions are addressed here 
for both thin and thick discs. Although quantities like velocity can always be formally 
separated into mean and fluctuating parts, the turbulence gives a negligible non-zero RMS 
error to the mean only when the disc radius is much larger than the scale of the turbulence. 
This RMS error will be estimated here as a precision measure of the a formalism, as a 
function of the averaging time. 

2. Precision of a-Accretion Disc Theory 

The usual slim disc equations are derived (e.g. Abramowicz et al. 1988; Narayan & Yi 
1995a) by vertically averaging the continuity, Navier- Stokes, and energy equations and with 
the magnetic field incorporated only as an additional pressure. Without presenting the 
formalism, here I assume the equations of e.g. Narayan & Yi (1995a) hold, but emphasize 
that the standard simple replacement of the micro-physical viscosity with a turbulent 
viscosity hides the requirement of radially and/or temporally averaging (in addition to the 
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usual azimuthal and vertical smoothing). For the radial average, a scale £ must be chosen 
such that l t b < £ < R, where R is the disc radius. The spatial average of a quantity like 
velocity V(i2) is then V (i2) = (V(R)) S ^ {V(R,t')}^ = Ja<£ V(il + A, £')<iA, where the 
similarity follows from the assumption that the time dependence is only due to turbulent 
fluctuations which are intended to be smoothed for mean quantities. The subscript 
indicates the mean quantity to be used in standard a disc theory. For a temporal average, 
taken over a duration t obs , we have V (R) = (V(R))t obs = (l/t obs ) j* +tobs V(R,t')dt' . 
The temporal average is meaningful only over times such that t Q b s > ttb where t t b is the 
dominant energy containing eddy turnover time scale. 

How precise is the assumption that the mean speed behaves as a steady monotonic 
function of R in the presence of turbulence? Note that "precision" is taken to mean 
that defined in Bevington (1969) with the a disc theory as the measuring device. A 
relative precision error (RPE) measures how effective the theory is at predicting, not how 
accurate the predictions are. The RPE error around the total mean speed Vq is given by 
AVq/Vq = [(AVb)/z/Vo + df>V (AR) rs /Vo], and the two terms on the right measure two 

1 /2 

RMS contributions to the RPE. The first can be approximated by (AVo)/z = vtb/NJ ~ 
Vtb{£,/hb + tobs/ttb) -1 ^ 2 where Nfi measures the "effective" number of eddies per radial 
averaging length. With increasing t b s , Nfi can well exceed the "snapshot" t D b s = 
value i/ltb- The second RMS contribution to the RPE above results from the fact that 

]/2 1/2 

R- (i/2)N rs LI " < R < R + (£ > /2)N rs L/ * is indistinguishable once £ is chosen. Here 
Nr}^ 2 ^ (1 + tobs/ttb) -1 ^ 2 and measures the "effective" number of averaging scales per £ 
which increases from its snapshot value of 1 for long t Q b s . The fluctuation and resolution 
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numbers (Nfi and N rs ) increase with t obs because the turbulence does not generate eddies 
in exactly the same location over time, so there is a smoothing effect which reduces the 
"effective" averaging and eddy scales. 

Since the total speed for both thin and thick discs is dominated by a contribution 
oc .R -1 / 2 , we can then estimate a total RPE for R as 

AR/R ~ [2(AV ) fl /V + (AR) rs /R] ~ [2(v tb /V )N- 1/2 + ^/2)N-^ 2 ]. (2) 

Assuming a constant accretion rate, (2) translates into an RPE in the luminosity given by 

AL V /(L V )~\4>\(AR/R) (3) 

where v is the frequency of emission and = |i?dfl[Ln(L„)]| which is usually 0.5 < < 
10 as will be addressed later. The RPE can be used to estimate the variability allowed for 
a given t obs . 

Though phenomenologically derived, the RPE formulae have properties which show 
that they capture the limiting cases correctly: First, for t obs » t t b, they are reduced as 
expected. Second, for t Q b s = 0, there is an optimal scale of 

i/R = Ut/R = (W^o) 2/3 (W#) V3 (4) 

for which the error is minimized: a larger £ reduces the RMS effect of the turbulent 
velocity, but one pays the price with a coarser spatial resolution. When ^ opt < kb, the 
RPE is dominated by the resolution term but then £ = l t b must be used since it is the 
minimum allowed. This will make the RPEs relatively independent of /?, since changing (3 
changes the relative importance of the two RPE terms while keeping (1) constant. 
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3. Energy Constraints and Relations Between Characteristic Speeds 

In a highly conducting turbulent plasma, the magnetic field is naturally amplified to 
the extent that v t b — B/(4irp) 1 ^ 2 = va, the Alfven speed (e.g. Parker 1979). Shearing 
box simulations, where turbulence is driven by a seed magnetic field (Stone et al. 1996; 
Brandenburg et al. 1995; Balbus & Hawley 1998), show va > v t b- Because of field line 
stretching, equipartition of turbulent and Alfven speeds is generally a more applicable rule 
of thumb than any relation between the particle and magnetic pressures. 

When the magnetic field is tangled on scales much smaller than those on which mean 
quantities vary (which for ADAFs likely requires temporal averaging, as seen below) aver- 
aging the Lorentz force gives an effective magnetic pressure 



Pmag/P ^ B 2 /(24np) = V \/6 = (1 - 0) P C 2 8 , 



(5) 



where (3 is a parameter. Using (5) in (1) and v t b — va we have 



l tb = aH/K{ 



(6) 



where K\ = 6(1 — (3). Because ltb < H, we have the constraint 



< 13 < (1 -« 2 /6). 



(7) 



The above relations also imply 



is 1 / 2 

V tb = A/ C s . 



(8) 



4. RPE of Thin Disc Models 
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For thin discs, H « R and V ~ V^, the Keplerian speed. From (3), (6) and (8) 
we have 

AL V /L V = \1>\AR/R ~ 2^|(WVb^)^ / ", 1/2 + 0.5|^| (e/^)^ 1/2 
„ 2|^|(g/fl)(^ 1 /l) 1 / 2 

" {Ma X [(l,22(^)V3(^)-2/3 ]+2 . 3xl 05 Us( ^ )( _a_ ) -l (] M ) -l ( _^_ ) -3/2 } l/2 

0.005|^|(a/0.01)(Ki/l)- 1 / 2 (g/i?)Max[l,22(X 1 /l) 2 / 3 (g/0.01)- 2 / 3 ] 
+ [l + 2.3xl0 5 U s (^i/l)(a/0.01)- 1 (M/M o )- 1 (i?/20i? s )- 3 / 2 ]V2 ' ^ 

where £ has been replaced by the Max[,] as per the discussion below (4), and relations 
(6), (8) and the definition of Nfi have been used. The temperature in an optically thick 
thin disc goes as T s oc R~ 3 / 4 (e.g. Frank, King & Raine 1992). Then, for example, in 
the Rayleigh- Jeans regime [hv « kT(R)) where L v oc u 2 J^ max T e (r)rdr, the luminosity 
within a radius R at a given frequency goes as L v oc z/ 2 i? 5 / 4 . Thus = 5/4. 

The RPE of (9) is small compared to what will be found for ADAFs. A careful check, 
keeping (7) in mind, ensures that for all allowed (3 the RPE < \i(j\H/R. (In section 7, we 
show that a and (3 can be related.) The RPE is further reduced for large t a b s . The low 
RPE results because H/R « 1 and v t b « v$ ~ vo^ for thin discs. 

5. Implications and RPE for Thick ADAF Discs 

For thick ADAF discs things are more subtle. From (7) and the standard ADAF 
choice of a = 0.3 (Narayan 1998a) we have < f3 < 0.985. Note when l t b = H the 
fact that H ~ R makes this upper limit extreme. Defining Ki = 2/(7 — 2(f) and using 
c s = Kl /2 V ff ~ (H/R)V ff (Narayan et al. 1998a), where V ff = (GM/R) 1 / 2 is the free-fall 

8 



speed, we have v t b = (KiK 2 ) 1 ^ 2 Vff, and thus from (6) and (8) 

hb/vtb = t tb = aH/(K lCs ) = aH/i^Ki^Vff = aR/^Vff). (10) 

Furthermore, defining K 3 = 12(1 — (3)/(7 — 2/3), the total mean speed is (Narayan 1998a) 
V ~ [(9a 2 /4)iff + Ks^^Vff, so we have 

K 4 = v tb /V = (K 1 K 2 ) 1 / 2 /[(9a 2 /4)^|+^3] 1/2 = [12(l-/3)(7-2/3)] 1 / 2 /[9a 2 +12(l-/3)(7-2/3)] 

(11) 

This is ~ 1 over the allowed range of < (3 < 0.985. Using (3) (4), (10) and (11) gives 

AL V /L V = ty\AR/R ~ 2ty\(v tb /V )NJ l 1/2 + \^/2)N-^ 2 

_ 1.22M(jr 4 /l)(<»/0-3) 1/2 

{Max[l,3.9(^)-§(^)§(^)^^)-i] + 2.4xl03^(^)(^)-i( 5 ^-)-i}-i 

0m\^\(a/0.3)(K 1 /3)-^(K2/0.58)^Ma4l,3.9(a/0.3)-i(K 4 /l)i(K 1 /3)^(K 2 /0.58)-^] 
[1 + 8 x 10 3 t o6s (a/0.3)-5(K 1 /0.3)(M/M G )- 1 ( J R/20 J R s )-i]-5 

(12) 

Recall that -^1,2,3,4 all depend only on (3. The rigorous form for ^ opt from (4) has also been 
employed. Over the allowed range < (3 < 1 — ct 2 /6, the RPE is relatively insensitive to 
(3 (though in section 7 we show that a and (3 may be related.) This is because decreasing 
(3 increases kb while lowering vtb and vice versa. The RPE is sensitive to both vtb and Itb- 

I now estimate \ip\ for various emission regimes based on ADAF scaling relations 
(Mahadevan 1997). Consider the radio Rayleigh- Jeans radio regime. Here L v ~ L„ c where 
v c is the peak frequency at each R < R max and is determined by synchrotron absorption 
In the ADAF, v c oc BT 2 oc T 2 R~ 5 / 4 , and is therefore a function of R. The spectrum traces 
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the envelope of peak frequencies, with each frequency corresponding to a particular R. For 
moderate accretion rates by ADAF standards (but below the critical value required for 
an ADAF solution) compressive electron heating is not important (Narayan et al. 1998a) 
and d R T e = 0. Using v = v c we have L v oc T^R~ 1 / 2 and = 1/2. When compressive 
heating is important, d R T e (R) ^ 0. Fitting the steadily falling T e (R) curve of Narayan et 
al., (1998b) I obtain LogT e (R) ~ 9.8 - 0.3Log(R/R s ) - 0m[Log(R/ R s )] 2 , so then 

Rd R [Ln T e (R)} = -0.7 - 0.3Log(R/ R s ), (13) 

and |V>| = 1/2 + 3.5 + 1.5Log(R/R s ). 

For v below v c = v c ^ min , the critical frequency corresponding to the maximum disc 
radius R m ax, the spectrum is simply oc v 2 at fixed R m ax, and in the constant T e regime, 
\ifj\ = 2. For the d R T e (R) ^ regime using (13), \ip\ = 1.3 + 0.3Log(R/R s ). 

In the Compton dominated sub-millimeter/X-ray regime L v oc v^ c T^ R~ 1 / 2 i>~ otc oc 
ji5+2a c ^-(2+5a c )/4^-a c ^ w ] 1 i c ] 1 j s sensitive to the Comptonization parameter a c and \ip\ = 

0.5 + 5a c /4 for constant T e {R). In the regime where (13) is applicable, \ip\ = 0.5 + 
5a c /4 + (5 + 2a c )[0.7 + 0.3Log(R/R s )]. In the Bremsstrahlung dominated sub-millimeter 
X-ray regime, L v oc -Ln(-R max /Rmir^FiT^T^Expl-hv/kT^ where F(T e ) is dominated 
by a term oc T e in the when kT e > m e c 2 and dominated by a term oc T^ 2 when kT e < 
m e c 2 . If d R T e = 0, \ip\ ~ l/Ln(R max /R min ) ~ 0.4, for R m ax/Rmin = 1000, but this 
is sensitive to radial dependences of T e since the latter appears in the exponential for 
the Bremsstrahlung regime. In the limit that T e > mc 2 , ~ 2.2(u/10 20 Hz)(T e /3 x 
lO 9 ^)- 1 [0.7 + 0.3Log(R/R s )]. For T e (R) < mc 2 , ~ | [0.5 - 2.2(v/10 20 Hz)(T e /3 x 
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10 9 K) - 1 ] [0.7 + 0.3Log(R/R s )] \ . 

Overall, for most regimes and reasonable (R/R s ), 1/2 < \tfj\ < 10. 

6. Implications for Observations of Presumed ADAFs 

We can see from (12) that for M ~ 1OM , (e.g. X-ray binary type systems) and 
\ip\ < 10, predictions probing the the inner 20R S and averages over t D&s > 10 3 sec are quite 
precise, that is AL U /L„ < 0.05. At R = 1000^, and |^| = 1/2, AL„/L U ~ 0.05. The 
allowed variability decreases as t~^/ 2 . 

Now consider the galactic nucleus of NGC 4258 with central mass M ~ 3.5 x 1O 7 M , 
For the radio emitting regime of this source near 22GHz, Herrnstein et al. (1998) found no 
detection of 22GHz emission in NGC 4258 with a 3-cr upper limit of 220/U Jy. This frequency 
is safely in the Rayleigh- Jeans regime and the best fitting models of ADAFs to NGC 4258 
have T e (R) approximately constant in this regime. Herrnstein et al. interpret this non- 
detection to mean that any ADAF proposed for this source (Lasota et al 1996) cannot 
extend outside radius defined by v c = 22GHz, namely R ~ 100R S . For the observations, 
tabs ~ 10 5 . But at 100i? s the t b s term does not contribute to significantly to (12). Since 
in this regime \ip\ ~ 1/2, we have from (12) AL^/L^ ~ 0.4, so this would reduce the 
significance of non-detection at 22 GHz to ~ la. 

For the Galactic centre, the presumed central mass is ~ 2.5 x 1O 6 M . Then from 
(12), at R = 20R S , t ofos must be > 10 3 seconds to contribute to significantly reducing 
the RPE. The X-ray observations and many radio observations above 10 GHz when taken 
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together provide enough total t obs for low RPE in this range (Narayan et al. 1998b and 
refs. therein) However for frequencies < 1GHz, R > 1000i? s , and the total t b s must be 
> 10 6 sec, for which there is insufficient data. Using \ip\ = 1.65 in (12), A.L V /L V > 1. 

Application of ADAFs to larger galactic nuclei (Fabian & Rees 1995) require much 
longer observation times for precise predictions. For M87, M ~ 3 x 1O 9 M so at 20i? s , the 
required t a b s time would be > 10 6 sec for the t obs term in (12) to reduce the RPE well below 
1, while at 1000i? s the the limit would be > 10 8 seconds of total time. X-ray observations 
have been made for 1.4 x 10 4 seconds (Reynolds et al., 1996) and radio observations have 
been made for only of order hours at particular frequencies, e.g. 2 x 10 4 sec at 1.7GHz 
(Reid et al., 1989). and 7.2 x 10 3 sec at 22 GHz (Spencer & Juror, 1986) More data are 
definitely needed in these sources. 

7. Conclusions 

The presumption that accretion discs are turbulent implies that the standard disc 
equations must be interpreted as mean field equations. Field line stretching equilibrates 
the turbulent and Alfven speeds so that the constraint on the ratio of magnetic to total 
pressures (Eq. 7) arises once it is assumed that the turbulent scale cannot exceed the 
disc height. For large a ADAFs this constraint is non-trivial and requires (3 < 0.985. For 
ADAFs, H/R ~ 1, so that 0.985 is really an extreme upper limit. 

The large turbulent scales and/or speeds for ADAFs lead to a relative precision error 
(RPE) or allowed variability in the predicted luminosity. Indeed, there are two contribu- 
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tions to the RPE. The first results from velocity fluctuations and the second comes the 
fact that the theory is unresolved on scales below the radial scale over which the averaging 
is taken. RPEs are larger ADAFs than for thin discs. The RPEs of (9) and (12) can 
be used to roughly predict allowed deviations from the a theory L v for a given t b s , and 
indicate when longer observations are needed to properly compare with disc models. The 
total RPE is reduced over large t a b s because such averaging amounts to smoothing over 
an ensemble of many turbulent realizations of the disc. 

The conclusions about the location of any transition radius to an ADAF in NGC 4258 
(Herrnstein et al. 1998) based on a 22GHz non-detection must be interpreted with this 
in mind since (12) suggests that an ADAF can significantly over-predict the luminosity. 
More observation time is then needed to fully rule out an ADAF outside 100i?s. For the 
Galactic centre below 1GHz, and for larger mass systems like M87 or M60 generally, more 
data than currently available are also needed to compare with ADAFs, since the RPEs 
only reduce significantly from 1 over long t b s - 

Finally, note that (9) and (12) do not necessarily presume that the Balbus-Hawley 
type (e.g. Balbus & Hawley 1998) shearing instability is the only source of viscosity. 
Making that assumption however, provides a 1-to-l link between a and (3: we have in the 
steady-state t tb ~ the instability growth rate, so t tb ~ R/vq^. Also, va ~ Vtb/3 1 / 2 from 
shearing simulations as discussed in section 3. Thus v = ac s H ~ v t bhb ~ v a(^/^ v o,<I>) 
which implies that a ~ 2(1-/3) for thin discs and a ~ 2(1-/3)(K 2 /^ 3 ) 1/2 = [2(1 -0)/3] 1/2 
for ADAFs. Eqns (9) and (12) are then sensitive to (3 mainly through the explicit a. 
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